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ABSTRACT 

We present the Spitzer/IRS mid-infrared spectrum of the infrared-luminous 
(L/zj = 4 X IO^^Lq) brightest cluster galaxy (BCG) in the X-ray-luminous clus- 
ter Z3146 {z = 0.29). The spectrum shows strong aromatic emission features, 
indicating that the dominant source of the infrared luminosity is star formation. 
The most striking feature of the spectrum, however, is the exceptionally strong 
molecular hydrogen (H2) emission lines, which seem to be shock-excited. The line 
luminosities and inferred warm H2 gas mass (~ IO^^Mq) are 6 times larger than 
those of NGC 6240, the most H2-luminous galaxy at z < 0.1. Together with the 
large amount of cold H2 detected previously (~ 10^^ Mq), this indicates that the 
Z3146 BCG contains disproportionately large amounts of both warm and cold 
H2 gas for its infrared luminosity, which may be related to the intracluster gas 
cooling process in the cluster core. 

Subject headings: galaxies: clusters: general — cooling flows — galaxies: cD — 
galaxies: active — infrared: galaxies 

1. Introduction 

Zwicky 3146 (Z3146)^ is an X-ray-luminous cluster of galaxies at a redshift of 0.29. Its 
soft (0.1-2.4 keV) X-ray luminosity ^ of 2 x 10^^ erg s~^ (Ebeling et al. 1998; Bohringer et 
al. 2000) is one of the largest known, and based on the ROSAT data, Z3146 was thought 
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to exhibit one of the most massive cluster coohng flows (Edge et al. 1994). Our analysis 
of Chandra data confirms that the gas cooling rate is indeed large with an estimated mass 
deposition rate of ~ 300 Mq yr~^ within a 50 kpc radius from the cluster center (Egami et 
al. 2006). 

The brightest cluster galaxy (BCG) is exceptionally active, showing a blue continuum 
and luminous nebular emission lines (Allen et al. 1992; Hicks & Mushotzky 2005) and an 
infrared luminosity of 4 x 10^^ Lq (Egami et al. 2006) . The mid-infrared spectral energy 
distribution (SED) suggests that the source of the large infrared luminosity is star formation 
at a rate of ~70 Mq yt^^, broadly consistent with the mass deposition rate mentioned above. 
This may indicate that the active star formation is caused by cooling cluster gas accreting 
onto the BCG, i.e., the cooling flow (Cowie & Binney 1977; Fabian & Nulsen 1977). At the 
same time, the BCG seems to contain a radio AGN that may contribute to its luminosity 
(Egami et al. 2006). 

To investigate the mid-infrared spectral properties in detail, we have observed the BCG 
in Z3146 using the Infrared Spectrograph (IRS; Houck et al. 2004) on board the Spitzer 
Space Telescope (Werner et al. 2004). 

2. Observations and Data Reduction 

The IRS spectra of the BCG in Z3146 covering 5-38 /xm were taken in the staring mode 
with the two low-resolution modules, Short-Low (SL; AA = 5.2 — 14.5/im) and Long-Low 
(LL; AA = 14.0 — 38.0//m), with a resolving power of ~ 64 — 128. The SL spectra were 
taken for 8 cycles of 60 s ramp durations for both the 2nd (SL2: 5.2-8.7 //m) and 1st (SLl: 
7.4-14.5 iim) orders while the LL spectra were taken for 8 cycles of 120 s ramp durations 
for the 2nd order (LL2: 14.0-21.3/im) and 6 cycles for the 1st order (LLl: 19.5- 38.0 /xm). 
The resultant integration times arc 960 s (SLl and SL2), 1920 s (LL2), and 1440 s (LLl), 
respectively, after coadding the spectra taken at the two nod positions. The slit width is 3'.'7 
for the SL and 10'.'7 for the LL module. 

Starting with the reduced and coadded two dimensional spectral image at each nod 
position, we performed sky subtraction in two passes, first by taking the difference between 
the two nod positions (spectral images targetting different orders were paired to avoid both 
the positive and negative spectra falling in the same order) and second by fitting and sub- 
tracting a linear fit to the background along each row in the cross-dispersed direction (only 
fitting the pixels well away from the source spectrum). Then, pixels with an abnormal signal 
near the source spectrum were manually masked out from the subsequent processing. From 
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the background-subtracted and cleaned spectral images, the source spectrum was extracted 
at each nod position using the program SPICE^ with extraction apertures reduced to 75% 
of the default widths. All the spectra were combined by resampling them with a uniform 
resolving power of 130 over the whole wavelength range. Finally, the combined spectrum 
was multiplied by 1.4 to match the IRAC 8.0 /^m and MIPS 24 /im photometry. 

3. Results 

The final combined IRS spectrum of the BCG in Z3146 is shown in Figure 1. The 
measured line fluxes and rest-frame equivalent widths are listed in Table 1. Along the 
spatial direction, the lines and continuum are unresolved (FWIIM~2 pixels), which sets an 
upper limit of 3'.'6 in SL and 10'.'2 in LL, respectively, on the size of the mid-infrared emitting 
region. This is consistent with the fact that the Z3146 BCG is barely resolved at 8 /^.m with 
a spatial resolution of 2" (Egami et al. 2006). 

Strong aromatic ("PAH") emission features arc detected at 6.2, 7.7, 8.6, and 11.3 fim. 
These features indicate that the mid-infrared (and therefore presumably the far-infrared) 
luminosity of this galaxy has a significant component of active star formation. 

The spectrum also shows exceptionally strong emission lines from molecular hydrogen 
(H2). We have securely detected the 0-0 S(l), S(2), S(3), S(4), and S(5) pure rotational 
lines. Although the 0-0 S(0) line is within the observed wavelength range {Xobs = 36.4/im), 
the noise in the LLl band increases dramatically beyond 35 /xm, preventing us from setting 
a meaningful upper limit on the line flux. The 0-0 S(6) line is blended with the bright 6.2 
/im PAH feature nearby. 

A few atomic fine-structure emission lines are also detected in the spectrum such as 
[Ne II] 12.8 iim and [Ne III] 15.5 /xm. The [S III] 18.7 /im hne is possibly detected as well, 
but because of the low signal-to-noise ratio, we treat it as a non-detection here. Also, no 
high-ionization hne (e.g., [Ne V] 15.6/xm) is seen. The small [Ne IIl]/[Ne II] line ratio of 0.2 is 
typical of a starburst galaxy (Thornley et al. 2000) although we cannot rule out a significant 
contribution to these neon lines from a shock-excited gas as discussed below. 

There are possible detections of the [Fe II] hnes at 5.3 and 26.0 //m. However, a longer 
integration time is needed to confirm the reality of these lines. 



^http://ssc. spitzer.caltech.edu/postbcd/spice. html 
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4. Discussion 

4.1. PAH Features 

In terms of the strength of the 6.2 //m PAH feature, the Z3146 BCG appears to be a 
typical starburst-dominated infrared luminous galaxy. The rest-frame equivalent width of the 
feature (0.67 fim) as well as its luminosity ratio to the total infrared luminosity (5 x 10"'^) are 
both within the range found for starburst galaxies, ~ 0.6 ±0.1 fira for the former (Weedman 
et al. 2005) and (6.3 ±3.2) x 10"^ for the latter (Peeters et al. 2004), respectively. Given the 
absence of any AGN-related features in the spectrum, the AGN contribution to the infrared 
luminosity is likely to be small. 

4.2. Molecular Hydrogen Lines 

The most striking feature of the spectrum is the exceptionally strong emission lines from 
molecular hydrogen. The luminosities of the pure rotational lines are on average 6.3 times 
larger than those of NGG 6240, the most H2-luminous galaxy at 2; < 0.1. 

The H2 line strengths are even more remarkable when compared with the infrared lu- 
minosity of this BCG. Figure 2a plots the line luminosity of the H2 0-0 S(l) emission line 
against the infrared luminosity of various galaxies. The H2 0-0 S(l) /infrared luminosity 
ratio of this BCG is extreme (0.25%), almost an order of magnitude larger than the ratio for 
NGC 6240 (0.03%), and a factor of 50 above the trend hue for the other galaxies (0.005%). 

Figure 3 shows the H2 level population in the Z3146 BCG. For comparison, the level 
population in NGC 6240 is overplotted after being scaled by a factor of 6.3. The level popu- 
lations are quite similar in shape between the two galaxies in the lowest energy states {v = 0; 
J — 3 — 7). These lowest energy states are almost always thermalized (i.e., coUisionally ex- 
cited) , and therefore the similarity suggests that although the total mass of warm H2 differs 
between the two galaxies, the fractional H2 mass distribution as a function of temperature 
is similar up to a gas temperature of ~960 K (i.e., the H2 excitation temperature derived for 
the V — 0, J — 6 and 7 states using the 0-0 S(4) and S(5) lines). 

Assuming that this similarity of the level population extends down to the v — 0, J — 2 
state, from which the 0-0 S(0) line originates, we estimate the total warm H2 mass in the 
Z3146 BCG to be ~ lO^^M©, i.e., 6.3 times larger than that of NGC 6240 (1.6 x lO^M© by 
Armus et al. (2006)). This mass is dominated by the ~160 K component, a temperature 
derived from the 0-0 S(0)/0-0 S(l) line ratio measured for NGC 6240. The cold H2 mass 
derived from the CO observation is also large, ~ IO^^Mq (Edge & Frayer 2003), but the 
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fraction of warm H2 mass is high, ~10%, similar to the 15% fraction found in NGC 6240 
(Armus et al. 2006). 

Figure 3 also shows that the H2 level population of the Z3146 BCG in the v — 1 
vibrational state is highly suppressed (by a factor of 30) with respect to that in the v — 
state when the NGC 6240 v — level population is used as a template. Although part of this 
suppression is due to the slit loss with the near-infrared spectroscopy as well as internal dust 
extinction, these effects, even when combined, are unlikely to be sufficient. For example, the 
size of the mid-infrared emitting region is 2" at most as already mentioned, so the I'.'S slit 
used to measure the f = 1 — near-infrared H2 lines (Edge et al. 2002) should capture most 
of the total line fluxes. Also, in this BCG, the extinction derived from the Balmer decrement 
is Ay — 0.6 mag (Crawford et al. 1999), and when corrected for this amount, the Ha-derived 
star formation rate agrees well with the infrared-derived star formation rate (Egami et al. 
2006), suggesting that the extinction cannot be much larger. In fact, near-/mid-infrared H2 
emission lines are known to show little extinction in infrared-luminous galaxies (Goldader et 
al. 1997; Higdon et al. 2006). 

This suppression of the v — 1 state suggests that the v — 1 level population in the Z3146 
BCG is subthermal as is the case with NGC 6240 (Egami 1998; Lutz et al. 2003). Although 
the cause of such a subthermal level population is not clear, one possibility is that the H2 
molecules are excited by a continuous (i.e., C-type) shock, which can produce subthermal 
level populations in the low vibrational states (e.g., Kaufman & Neufeld 1996; Le Bourlot 
et al. 2002). Note that alternative excitation mechanisms such as a jump (J-type) shock or 
far-UV radiation would produce a. v = 1 level population at or above the thermal level with 
respect to the v — state (e.g.. Burton 1992). Figure 3 shows that only a factor of ~7 scahng 
is required to bring up the v — 1 level population in the Z3146 BCG to that of NGC 6240. 
This correction factor is much easier to account for than the factor of 30 mentioned above, 
which would be required if the v — 1 state is thermalized. Shock excitation is also consistent 
with the previous near-infrared spectroscopic studies of nearby cooling-flow cluster BCGs 
(e.g., Jaffe et al. 2001; Wilman et al. 2002). 



4.3. Atomic Fine-Structure Lines 

The [Ne II] 12.8 pim line of this BCG is also strong. Figure 2b plots the [Ne II] 12.8 //m 

himinosity against the infrared luminosity for a sample of nearby infrared-bright galaxies. 
As shown by Sturm et al. (2002), there is a good correlation between the two luminosities 
among the starburst and Seyfert galaxies, and the Z3146 BCG again stands out by having 
an exceptionally strong [Ne II] 12.8 /^m line for its infrared luminosity (0.75%), 
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This overluminous [Ne II] emission line may also arise from a shock. However, the C- 
type shock invoked above to explain the H2 level population cannot produce a substantial 
amount of the [Ne II] line because the shocked gas in a C-type shock is largely neutral and 
molecular. To produce the [Ne II] line, the shock must be a fast (> 50 km s~^) J-type shock, 
that can dissociate and ionize the shocked gas (e.g., Hollenbach & McKee 1989). The Z3146 
BCG has an Ha velocity width of ~ 600 km s~^ (Allen et al. 1992), so such a fast shock is 
plausible. 

This probably means that various types of shocks may coexist in the Z3146 BCG. A 
similar conclusion was reached for NGC 6240 by Lutz et al. (2003) to account for the strong 
[O I] 63 iim line, which is also a sign of a J-type shock. Indeed, such a coexistence of the 
C- and J-type shocks is also seen in some Galactic outflow sources, in which both strong H2 
and [Ne II] lines are detected (e.g., Lefloch et al. 2003). Shock models can produce a wide 
range of [Ne III]/[Ne II] hne ratios depending on various shock conditions, and therefore 
it is possible to produce the starburst like [Ne IIl]/[Ne II] ratio with a strong shock (e.g., 
Binette et al. 1985; Hartigan et al. 1987). [Fe II] hues may also be used to constrain the 
shock conditions if their strengths are measured accurately. 

With respect to the [Ne II] 12.8 //m line, the [S III] 18.7 iim line is abnormally weak 
([S IIl]/[Ne II] <0.05). Such a low [S IIl]/[Ne II] ratio is often seen with ultraluminous 
infrared galaxies (ULlRGs), indicating that the gas phase abundance of sulfur is significantly 
below solar (Genzel et al. 1998). Starburst galaxies also show this sulphur "deficiency" 
although the [S IIl]/[Ne II] ratios are less extreme, and this may suggest that the sulphur is 
depleted onto dust grains in high-metallicity star- forming systems (Verma et al. 2003). 

5. Conclusions 

We present the Spitzer/IKS mid-infrared spectrum of the infrared-luminous BCG in 
the X-ray luminous cluster Z3146. The strong aromatic features indicate that the dominant 
source of the infrared luminosity is star formation. The most striking feature of the spectrum, 
however, is the exceptionally strong H2 emission lines. To our knowledge, this BCG has the 
most luminous H2 pure-rotational emission lines and most massive warm H2 gas component 
(~ lO^'^ Mq) known. Together with the large amount of cold H2 detected previously (~ 
IO^^Mq), this indicates that the Z3146 BCG contains disproportionately large amounts of 
both warm and cold H2 gas for its infrared luminosity. Since the parent cluster Z3146 has an 
exceptionally large gas cooling rate in the core, these massive warm and cold H2 components 
may be related to the intracluster gas cooling process. 
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Tablc 1. Z3146 BCG line fluxes 



Line/Feature 


-^rcst 


Flux 






(urn) 


(10-1^ W m-2) 


(/xm) 


PAH 


6.2 


3.0±0.2 


0.67±0.05 


H2 0-0 S(5) 


6.9 


2.1±0.1 


0.49±0.02 


[Ar 11]'^ 


6.99 


~0.6 


-0.15 


PAH 


7.7 


6.9±0.2 


1.72±0.06 


H2 0-0 S(4) 


8.0 


0.8±0.1 


0.19±0.02 


PAH 


8.6 


1.2±0.2 


0.31±0.04 


H2 0-0 S(3) 


9.7 


2.3±0.1 


1.72±0.05 


PAH 


11.3 


1.3±0.1 


0.47±0.07 


H2 0-0 S(2) 


12.3 


1.2±0.1 


0.26±0.02 


[Ne 11] 


12.8 


4.2±0.1 


1.34±0.03 


[Ne 111] 


15.6 


0.8±0.1 


0.30±0.04 


H2 0-0 S(l) 


17.0 


1.4±0.1 


0.55±0.04 


[S 111] 


18.7 


< 0.4 


< 0.2 



'^The [Ar 11] 6.99 /im line was not detected directly, 
but its presence was set to the maximum level consistent 
with the S(5) line shape. 
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Fig. 1. — Combined IRS spectrum of the BCG in Z3146. The two points with error bars 
show the IRAC 8.0/i.m and MIPS 24 /xm photometry from Egami et al. (2006). The IRS 
spectrum has been scaled (xl.4) to match these photometry points when the spectrum is 
integrated over the corresponding passbands. 



- 12 - 




LiR (L©) LiR (Lg) 

Fig. 2. — (a) H2 0-0 S(l) 17 /iin line luminosity vs. infrared luminosity (Lj^). The solid and 
dashed lines correspond to i^o-os(i) = 0.005 and 0.25% of Lj^. The H2 line luminosities were 
taken from Armus et al. (2006) (NGC 6240), Rigopoulou et al. (2002), and Higdon et al. 
(2006). Note that the H2 luminosity correlates with the infrared luminosity, so the sample 
studied by Higdon et al. (2006), which extends beyond a redshift of 0.1 and therefore includes 
more infrared-luminous galaxies, contains galaxies more H2 luminous than NGC 6240; (b) 
[Ne 11] 12.8 iim line luminosity vs. Lir. The sohd and dashed lines correspond to -L[jve7/] = 
0.06 and 0.75% of Ljji. The [Ne II] line luminosities were taken from Genzel et al. (1998) 
while the infrared luminosities were taken from Sanders et al. (2003). 
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5.0»10^ I.O'IO^ 1.5»10^ 

EA(K) 

Fig. 3. — H2 excitation diagram showing the level population of warm H2 in the Z3146 BCG. 
The level population in the ground vibrational state {v — 0) is based on our observation while 
the V — 1 level population was derived from the data by Edge et al. (2002). The X axis 
is the H2 energy level expressed in temperature {Ei is the energy level of the i-th state; k 
is the Boltzmann constant). In the Y axis, Ni refers to the total number of H2 in the i-th 
energy state while gi is the statistical weight of that state (ortho-para ratio was set to 3:1). 
Assuming that the H2 lines are optically thin, Ni was calculated as Ni = Li/lAihui), where 
Lj is the observed line luminosity, Ai is the Einstein coefficient for that transition, h is the 
Planck constant, and Ui is the frequency of the line. In comparison, the H2 level population 
of NGC 6240 is also shown based on the published data (van der Werf 1996; Sugai et al. 
1997; Lutz et al. 2003; Armus et al. 2006). The data from Sugai et al. (1997) were plotted 
by assuming that the 1-0 S(l) fine flux agrees with that of van der Werf (1996). No other 
scaling was applied. The solid line connects the v = and 1 states in the Z3146 BCG while 
the dotted line connects the v — states in NGC 6240. 



